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Abstract - We hnd excitations lower in energy than known phonon modes in underdoped 
La 2 -ccSrccCu 04+5 (x=0.08), with both inelastic X-Ray scattering (IXS) and inelastic neutron 
scattering (INS). A non dispersive excitation at 9 meV is identihed and is also seen by INS in 
(La,Nd) 2 -ccSra;Cu 04 + 5 , with 40% Nd substitution. INS also identifies a still lower energy disper¬ 
sive mode at low q in the Nd free sample. These modes are clearly distinct from the longitudinal 
acoustic phonon and correspond in energy to the Zone Centre modes measured by optical spec¬ 
troscopy and associated with stripe dynamics. 


The relevance of exotic charge-modulation states and their co-existence with the super¬ 
conducting state in underdoped high temperature superconductors is a subject of current 
debate. This state has now been conclusively detected in a variety of materials where this 
ordering is statie and is linked to lattice deformations and strain Bfl Reports of static 
stripes exist also in materials like nickelates |8 -10 where superconductivity is never attained 
while in cuprates pinning of stripes by lattice deformations is accompanied by suppression 
of superconductivity (2]-[m. Available evidence of dynamical charge-modulations coexisting 
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with superconductivity is mostly indirect |llf[l4| . 

A related debate concerns the interpretation of charge-modulations in cuprates as either 
stripes, z.e., driven by Coulomb repulsion [^, or as Charge-Density-Waves (CDW) either 
driven by lattice deformations as originally claimed 14 , or, more recently, by a nesting of the 
Fermi Surface 


7[|15| . The two interpretations give completely different doping dependences 
of the charge modulation wave-vector: for the first one we expect an increase with hole¬ 
doping [2||^ , while for the CDW nesting-driven case we expect a decrease with hole-doping 
15 . Both scenarios yield comparable wave-vectors close to x 1/8 where they are most 
easily detected 16 . Hence it is fundamental to characterize possible signatures of such 


charge-modulations across the entire phase diagram of cuprates. 

Optical spectroscopies such as Raman and infrared experiments have revealed the pres¬ 
ence of two new low energy modes in underdoped La 2 -a;Sra,Cu 04+<5 ^T\ . The first mode 

occurs at 25 cm“^ 3 meV) at low temperature (iTpSpol and shows a strong temperature 
dependence: as the temperature increases the mode frequency increases up to ~ 200 cm“^ 
while its intensity decreases, eventually becoming zero above ~ 200 K 17 . The second 
mode occurs at higher energy, at about 80 cm“^ 10 meV) at low temperature. More 
recently, time-resolved optical spectroscopy in the femto-second regime showed a resonant 
mode of about 2.0 THz frequency that corresponds to the same energy range 70 cm“^ or 
8 meV) [^. Such mode is interpreted as an amplitudon, in line with Ref. [^, suggesting 
a charge-density-wave like model. 

The interpretation of these Raman modes is the subject of debate: do they result solely 
from scattering by electronic excitations 17 19 or are they coupled charge-lattice modes 


18 


Recently IXS experiments have reported anomalies in low energy acoustic modes or 
additional modes at lower energies in cuprate Bi2201 23 and oxychlorides 24 . In the 


25 


”123” family (as YBCO, NdBCO,...) acoustic phonon anomalies have been detected 
and linked to CDW [7 25 . However the charge modulations in the ”123” family bear 
distinct feature from the one of the ”214” family {e.g. (La, Nd) 2 -xSrxCu 04 ) we study here. 
In particular the ”123” family show no associated spin-fluctualtion in strike contrast to 
what observed in the ”214” one [^[^. 

Here we report an observation of anomalous low energy modes for the first time in un¬ 
derdoped La 2 -a:Sra,Cu 04+<5 using both IXS as in the above cited studies, as well as INS, 
thus narrowing drastically the possible origins of the nature of the excitation and its disper¬ 
sion. We identify two additional modes at energies compatibles with the results obtained 
by optical spectroscopies 17 18 20 at the zone centre, and extend the observation to the 


whole Brillouin Zone, as well as to the Nd substituted sample. In the y=0.4 Nd substituted 
sample where the charge modulation is pinned we still observe a non dispersive mode 
while the lower energy dispersive mode is much weaker, and even not detected for reduced 
wave-vectors (g, 0, 0) with 0.2 < q < 0.4. Using inelastic neutron scattering, we have ex¬ 
plored the same (Q,E) region in a crystal with the same doping, x=0.08, in order to identify 
the possible origin of these excitations by comparing neutron and non-resonant X ray cross 
sections. We find strong evidence that they are coupled to lattice modes. 

IXS and INS measurements were performed on two underdoped La 2 -a:Sra,Cu 04 +^ both 
with x=0.08, where Raman measurements report the strongest intensities for the extra 
modes. The crystal for INS measurements (sample label SI) has cylindrical shape with an 
average diameter of about 6 mm, a height in the order of of ^ 35 mm and Tc=20 K, while 
on IXS experiments we measured a second crystal with a platelet shape, about 0.5 mm thick 
and few mm wide, with Tc=(20.7±0.2) K (sample label S2). With INS we also measured a 
La 2 -a:-yNd^Sra,Cu 04 +^ crystal of cylindrical shape with an average diameter of about (5 ± 
0.5) mm and a height of ^ 25 mm grown from the melt in an image furnace by the traveling 
solvent floating zone method, with with x=0.08 (as above) and y=0.4 (Tc=(7.5±0.25) K, 
sample label SNd). 

The sample was mounted in reflection geometry on the cold finger of a closed-loop helium 
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cryostat. The sample c crystal axis was perpendicular to the scattering plane, with a beam 
spot ofh X v~0.3 X 0.1 mm, while the cm^ sample for the INS experiment was completely 
irrigated by the neutron beam. The highly collimated X-ray beam and small sample size 
result in a small probed volume and an angular width of 0.01° FWHM for the (4 0 0) 
Bragg reflection. In the neutron experiment we measured 1.5° (SI) and 0.8° (SNd) FWHM 
for the (2 0 0) reflection. Phonons were measured in the fourth Brillouin zone in the 
IXS experiment and from second to third zones by INS. The IXS spectrometer 
was used with simultaneous measurements from 5 analyzers, and a resolution of 
(see e.g 


Ref. for general details, and Ref. 16 


3meV 

for the setup specific to the present 
experiment). We note that IXS spectra do not have contributions from incoherent scattering 
or multiple scattering to their S{Q^uj) dynamical structure factor. For X-rays incoherent 
Compton scattering occurs above ^ 1 eV energy transfer due to the small electron mass, 
while multiple scattering is strongly diminished by the high photoelectric absorption cross 
section of cuprates. The multiphonon part is then typically the only intrinsic background 
in IXS for meV energy transfers, and can be easily singled out 29 . On INS experiments 


we occasionally observe additional contributions in the measured energy range, mostly from 
higher monochromator harmonics. These can be calculated and thus easily eliminated. 
A standard lattice dynamical calculation using a shell model code 30 was performed to 
calculate the phonon dispersion and the expected intensities from inelastic scattering as 
described in Ref. 


16 


In the left panel of Fig. [^we plot the IXS intensity as a function of energy for reduced 
wave-vectors q = (0.387,0.045,0) and q = (0.4, 0,0). We detect three phonon modes which 
can be identified as the longitudinal acoustic mode (labeled LA, 16 meR), the first optical 
mode (labeled ILO, ^ 20 meV) and the second optical mode (2LO, ^ 35 meV). These are 
expected from shell model calculations but also from earlier measurements [m|[^ . We also 
observe an additional low energy mode (labelled E, 9 meV) which at this wave-vector 
point clearly emerges as a peak with intensity comparable to the other phonon modes. The 
deviation from the longitudinal condition along (100) for the q = (0.387,0.045,0) spectra is 
negligible according to calculation, with less than 3% transverse contribution. Indeed, the 
data shows no visible difference with the spectra measured in exact longitudinal condition 
with q = (0.4, 0,0). We also check for contribution coming from other symmetries such 
as transverse modes. Though the measuring geometry used forbids the measurement of 
transverse phonons, the non-zero solid angle covered by the analyzer may favor such a 
contribution. However in our case the analyzer acceptance angle is too small to account for 
such a parasitic signal from transverse phonons. In Fig. left panel, we show the calculated 
IXS intensity, integrated over the whole scattering volume in reciprocal space, according to 
our resolution and mosaic spread, and convoluted with the experimental energy resolution, 
by the blue filled area The calculated intensity is concentrated under the measured 
phonon modes, but there are no contribution from phonon modes to the scattered intensity 
^ 9meV. The right panel of Fig. shows representative IXS and INS data at g = (0.4, 0., 0) 
for temperatures ranging from 24 to 300 K, together with the corresponding fits using a 
model function as described elsewhere 24 34 . INS data are fitted to a damped harmonic 


oscillator convoluted with the instrumental function as calculated by a program including 
the spectrometer characteristics 35 . For all temperatures, we detect the phonon modes 
and the additional low energy excitation labelled ”E” in Eig[^ left panel. A second mode, 
seen by Raman at lower energies, is also observed in the INS data, but only in a limited 
range of the Brillouin Zone, from q = (0.07,0,0), where it intercepts the acoustic phonon 
mode dispersion, at approximately 5 meV, to g = (0.4, 0,0), where it seems to merge with 
the higher mode ”E” at about 10 meV. 

Both for IXS and INS we do not have low energy data close to the Bragg reflection. 
However, we have zone centre data around forbidden reflections {2H + I, 0, 0). We represent 


similar procedure has been used in Ref. l33|, and described in the supplemental material, Fig. FigSI2. 
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Q=(3,0,0)+q; T=24K 



Q=(3 0 0)+/-(0.4,0,0) 



Fig. 1: (Color online). Left panel. Energy loss IXS spectra in longitudinal configuration for 
Q = (3 + ^,^>,0) at T = 24 K. Black circles: ^ = 0.387 and v = 0.045; red (gray) circles: ^ = 0.400 
and V = 0.0. The blue (gray) filled area shows a calculation of the expected intensity from our model, 
integrated over the q-resolution volume around Q = (3.4,0,0), convoluted to the instrumental 
energy resolution function. Right panel. Energy loss IXS and INS spectra and best fit to a sum of 
harmonic oscillators convoluted to the instrumental energy resolution function. The data are taken 
at Q = (3.4, 0, 0), for IXS, as in the left panel, and at the equivalent position Q = (2.6, 0, 0) for INS 
data. The latter are taken for different temperatures and sample of La 2 -cc-yNdySra,Cu 04+5 with 
x=0.08 and different Nd substitution: y=0 for sample SI and S2 and y=0.4 for SNd. 


our data in an extended Brillouin Zone, where point with reduced wav-vector (g, 0, 0) for 
q > 0.5 are unfolded up to these forbidden reflections, as usually done for I4/mmm-like 
symmetries 36 . But these points in fact can be folded back to the Zone Centre, and have 


been measured in the INS experiment, as shown in Fig. and can be directly compared to 
the optical results. In that figure we show indeed how the energy found in the Raman results 
of Ref. 18 directly compare to our INS data, giving strong evidence that we are actually 


measuring the same feature. Unfortunately, such measurements are not possible with IXS, 
because of the presence of a strong parasitic elastic scattering, while the INS measurements 
are limited to E > 5 meV, due to geometrical constraint at large exchanged wave-vector, 
so that we can access only the excitation E, identified as an amplitudon mode in e.g. Ref. 

18 22 , while the feature at low energy can not be measured at the zone centre. 


In Fig. [^we compare INS intensities for several fixed-energy g-scans (close to the Bril¬ 
louin Zone Centre) and fixed-Q energy scans on the same sample at low temperature (10 
K) with phonon modes expected from calculations. For energy losses below 5meV in the 
g-scans of the left panel only the acoustic phonon is seen. The energy scans of the right 
panel clearly show a strong anomalous mode at ~ 9meV. 

In Fig. [^we show the experimental and calculated dispersion of phonons in the extended 
zone in the usual representation for these systems 36 . The experimental dispersion is 


obtained from that of the fitted peaks in the experimental energy loss spectra as shown in 
Fig. [^and|^ Phonon energies and their dispersion correspond well to the calculation in the 
second half of the extended Brillouin zone, from (tt/u, 0 0) to (27r/a, 0 0) while the extra 
mode already shown in Fig. and now shows up as a fiat, optical-like band. In the first 
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Fig. 2: (Color online). Raw INS energy scan on La 2 -xSrxCu 04+5 with x=0.08 (SI) at Q = 
(3,0,0) and T=150 K. Red lines shows a best fit to a sum of harmonic oscillators convoluted 
to the instrumental energy resolution function, black lines representing each harmonic oscillators 
contribution (A single black line when only one contribution is fitted). Blue vertical line indicates 
the calculated position for the acoustic (low energy) and first optic (high energy) mode, while the 
black, dashed line indicate the position of the Amplitudon” according to Raman data in Ref. p^ . 


half of the Brillouin zone this mode interacts strongly with the acoustic phonon opening a 
gap as in an avoided crossing (or anti-crossing) (see middle figure of the right panel of Fig. 
[^. A lower energy dispersive component also appears which increases in energy with g, 
finally merging into the higher energy anomalous mode at g = (0.4,0,0). 

The INS measurements confirm the observations made using IXS and we clearly detect 
the additional mode at approximately 9-10 meV. Moreover we observe its interaction with 
the acoustical phonon as well as another anomalous dispersive mode at lower energy. The 
IXS and INS frequencies superpose well in the reduced q-range where we have access to 
both. This additional mode is consistently seen with different experimental set-ups and 
techniques, over different Brillouin zones, and on samples with different mosaic. This has 
several important implications: 


a) Any possible instrumental artefact is excluded; 

b) The anomalous modes are linked to the periodicity of the crystal since measurements 
were made in different Brillouin zones; 


c) The mode ”E” implies a displacement of the centre of the mass of the atoms, as it is 
seen both in the INS and IXS channel. In any case, a pure charge (in principle visible 
only by IXS) or spin (visible only by INS) excitation is excluded. 


In fact, if this peak was of purely electronic origin we would not have been able to detect 
it with INS, and even for IXS the corresponding cross section would be several orders of 
magnitude smaller than that of the low energy phonon modes. 

We also confirm that the additional modes appear only along (1 0 0), but not along the 
diagonal (110) nor in transverse configuration, where the data match previously reported 
experiments. In addition, we observe a clear anti-crossing with the acoustic longitudinal 
mode of A1 character, indicating that the LA and the ”E” mode possess the same symmetry. 

Erom a phenomenological point of view, parallels can be drawn between the symmetry 
breaking associated with stripe formation and that coming from the formation of charge 
density waves. This point has been recently studied by optical, time-resolved experiment 
22 in the same system, as well as in a very similar experiment in YBCO 25 . In the latter 
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q (r.I.u.) 


INS scans SNd Q=(2 0 0)+q T=10 K 



Fig. 3: (Color online). Raw INS data of La 2 -x-yNd^SrccCu 04+5 with x=0.08 and y=0.4 (SNd) along 
Q = (2 + g, 0, 0) at ~ 10 K. Left panel: fixed-energy, g-scans, close to the Zone Centre (2, 0, 0). Right 
panel: fixed-g energy scans. Red lines shows the best fit to a sum of harmonic oscillators convoluted 
with the instrumental energy resolution function; black lines represent individual harmonic oscillator 
contributions (A single black line when only one contribution is fitted). Blue vertical dashes indicate 
the calculated position for the acoustic mode at ±q in the left panel, and acoustic (low energy) and 
first optic (high energy) mode in the right panel. 


case, however, a true static lattice deformation accompanies the electronic transition and has 
a drastic effect on the phonon spectrum: softening of phonon branches to the point where 
the static lattice symmetry is altered. This softening and the new symmetry generate two 
new modes: the acoustic like phason with a very rapid dispersion given by the softening, and 
the optical like amplitudon 37 . In our case we do not observe similar softening, associated 
to a possible ”phason” mode, even when we add Nd, that is supposed to pin the ’’stripe” 
order. On the one hand the effects may be subtle since fluctuating stripes would not bring 
about a change in the static lattice symmetry. On the other at a Sr-doping x quite far from 
1/8, the pinning effect could be weak. 


Our measurements point to a local symmetry breaking on the time-scale corresponding 
to THz frequencies, resulting from the coupling of the lattice with spatial charge modula¬ 
tions, possibly with dynamical and short range character. They correspond favorably with 
the Brillouin Zone Centre measurements by Raman and infrared experiments for a similar 
doping, where they are identified as due to stripes. An important open question has been 
the possible link between this phase and the lattice. We show here that the low energy 
lattice dynamics bears an observable fingerprint of this phase. Further studies close to the 
X « 1/8 doping where the Nd substitution pins the charge-modulations as well as in the 
tetragonal phase far from structural transition will allow us to disentangle possible lattice 
effects from purely electronic ones. 
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Fig. 4: (Color Online). Dispersion of the longitudinal phonon modes with A symmetry along 
(3+^ 0,0). Note that this is an extended Brillouin Zone, with the zone centre at (2H,0,0) and zone 
boundary at (2H+1, 0, 0) corresponding to the in-plane (27r/a , 0) point. Calculations (lines) are 
compared to the experimental results, as extracted from IXS (red filled diamonds) and INS (blue 
and white filled circles) spectra. Data are shown up to the 2LO mode, the complete set of optic 
modes is reported in Ref. [^. INS data are reported for two different Nd substitutions as in Fig. 

Raman data are taken from Ref. [^, and we trace a dashed line connecting these data at the 
Brillouin Zone Centre with the same ones unfolded at the extended zone boundary. 
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